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Abstract 

As generalizations of results of Christandl et al.[l, 2] and Facer et al.[3], Bernasconi 
et al.[4, 5] studied perfect state transfer (PST) between two particles in quantum net- 
works modeled by a large class of cubelike graphs (e.g., the hypercube) which are the 
Cayley graphs of the elementary abelian group Z%- In Refs. [6, 7], respectively, PST of 
a qubit over distance regular spin networks and optimal state transfer (ST) of a d-level 
quantum state (qudit) over pseudo distance regular networks were discussed, where the 
networks considered there were not in general related with a certain finite group. In this 
paper, PST of a qudit over antipodes of more general networks called underlying net- 
works of association schemes, is investigated. In particular, we consider the underlying 
networks of group association schemes in order to employ the group properties (such as 
irreducible characters) and use the algebraic structure of these networks (such as Bose- 
Mesner algebra) in order to give an explicit analytical formula for coupling constants in 
the Hamiltonians so that the state of a particular qudit initially encoded on one site will 
perfectly evolve to the opposite site without any dynamical control. It is shown that the 
only necessary condition in order to PST over these networks be achieved is that the 
centers of the corresponding groups be non-trivial. Therefore, PST over the underlying 
networks of the group association schemes over all the groups with non-trivial centers 
such as the abelian groups, the dihedral group Dm with even n, the Clifford group CL(n) 
and all of the p-groups can be achieved. 

Keywords: Perfect state transfer, Qudit, Underlying graphs, Group asso- 
ciation schemes, Stratification 

PACs Index: 01.55.+b, 02.10.Yn 
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The quantum communication between two parts of a physical unit, e.g., a qubit, is a crucial 
ingredient for many quantum information processing protocols [8] . There are various physical 
systems that can serve as quantum channels, one of them being a quantum spin system. In 
view of applications like the communication between registers in quantum devices, the study 
of natural evolution of permanently coupled spin networks has become increasingly important. 
A special case of interest consists of homogenous networks of particles coupled by constant and 
fixed (nearest-neighbor) interactions. An important feature of these networks is the possibility 
of faithfully transferring a qubit between specific particles without tuning the couplings or 
altering the network topology. This phenomenon is usually called perfect state transfer (PST). 
Quantum communication over short distances through a spin chain, in which adjacent qubits 
are coupled by equal strength has been studied in detail, and an expression for the fidelity 
of quantum state transfer has been obtained [9, 10]. Similarly, in Ref. [11], near perfect 
state transfer was achieved for uniform couplings provided a spatially varying magnetic field 
was introduced. After the work of Bose [9], in which the potentialities of the so-called spin 
chains have been shown, several strategies were proposed to increase the transmission fidelity 
[12] and even to achieve, under appropriate conditions, perfect state transfer [1, 2, 13, 14, 
15, 16, 3]. Recently, A. Bernasconi, et al. [4] have studied PST between two particles in 
quantum networks modeled by a large class of cubelike graphs. Since quantum networks (and 
communication networks in general) are naturally associated with undirected graphs, there is a 
growing amount of literature on the relation between graph-theoretic properties and properties 
that allow PST [5]. In Ref. [6], the so called distance-regular graphs have been considered 
as spin networks (in the sense that with each vertex of a distance-regular graph a qubit or a 
spin 1/2 particle was associated) and PST over them has been investigated. In Ref. [17], state 
transfer over spin chains of arbitrary spin has been discussed so that an arbitrary unknown 
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qudit be transferred through a chain with rather a good fidelity by the natural dynamics of 
the chain. In the recent paper [7], the authors have investigated optimal state transfer (ST) 
of a (i-level quantum state (qudit) over pseudo distance regular networks, where it was shown 
that only for pseudo distance regular networks with some certain symmetry (mirror symmetry) 
in the corresponding intersection numbers (consequently their QD parameters), PST between 
antipodes of the networks can be achieved. 

In the present paper we will consider the more general graphs that are the underlying 
graphs of group association schemes [18] and give necessary and sufficient conditions for PST 
of a qudit state in quantum networks modeled by a large class of graphs with group structure. 
An association scheme is a combinatorial object with useful algebraic properties (see [19] for 
an accessible introduction). The theory of association schemes has its origin in the design of 
statistical experiments. This mathematical object has very useful algebraic properties which 
enables one to employ them in algorithmic applications such as the shifted quadratic character 
problem [20] (in this problem, the association scheme is the Paley scheme which corresponds to 
a strongly regular graph, the Paley graph). A d-class symmetric association scheme (d is called 
the diameter of the scheme) has d + 1 symmetric relations Ri which satisfy some particular 
conditions. Each non-diagonal relation Ri can be thought of as the network (V,Ri), where 
we will refer to it as the underlying graph of the association scheme (V is the vertex set of 
the association scheme which is considered as vertex set of the underlying graph). In fact, an 
association scheme partitions the relationships between pairs of vertices into classes, so that for 
an arbitrary chosen vertex (referred to as reference vertex), one can stratify the vertices into 
distinct classes or stratas according to its relationships with all of the other vertices. Moreover, 
this stratification is independent of the choice of reference vertex. In the problem of transfer 
of an arbitrary qubit state which is considered in this work, we are given with iV spin-1/2 
particles as the corresponding vertex set; Then, for a given particle associated with a vertex of 
the underlying graph, the strength of its interaction with other N — 1 particles is determined 
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according to its relationship with the other vertices defined via the relations of the correspond- 
ing association scheme; In Refs.[21], [22], algebraic properties of association schemes have been 
employed in order to evaluate the effective resistances in finite resistor networks, where the 
relations of the corresponding schemes define the kinds of resistances or conductances between 
any two nodes of the networks. In Ref.[23], a dynamical system with d different couplings 
has been investigated in which the relationships between the dynamical elements (couplings) 
are given by the relations between the vertexes according to the corresponding association 
schemes. Indeed, according to the relations the so-called adjacency matrices Ai are defined 
which form a commutative algebra known as Bose-Mesner (BM) algebra. One of the important 
preferences of association schemes is their useful algebraic structures that enables one to find 
the spectrum of the adjacency matrices relatively easy; Then, for different physical purposes, 
one can define particular spin Hamiltonians which can be written in terms of the adjacency 
matrices of an association scheme so that the corresponding spectrums can be determined 
easily. Group association schemes are particular schemes in which the vertices belong to a 
finite group and the relations are defined based on the conjugacy classes of the corresponding 
group. Working with these schemes is relatively easy, since almost all of the needed informa- 
tion about the scheme, for instance the so-called eigenvalue matrices P and Q associated with 
the scheme, can be obtained via the character tables of the corresponding groups. We will 
use the technique of the stratification of the underlying networks of group association schemes 
(this technique can be used even for some graphs which are not the underlying graphs of as- 
sociation schemes [24, 25], not only for the purpose of state transfer, but also in investigating 
the continuous time quantum walk (CTQW) over the undirected graphs [26]- [30]) and employ 
their algebraic structures in order to calculate the transition probability amplitude between an 
arbitrary chosen reference vertex (the identity element of the corresponding group is associated 
with this vertex) and the antipode vertex (any element of the center of the group which forms 
one element conjugacy class or one element strata, can be considered as the corresponding 
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antipode vertex) and optimize it in order to attain PST between them (transmission fidelity 
attains to maximum value one). As we will see, the preference of this employment is that we 
are able to give analytical formulas for coupling strengths in particular Hamiltonians imposed 
to these graphs, in terms of the irreducible characters of the corresponding group, in order 
for PST to be achieved. In fact, we show that for those such networks in which the group G 
has non-trivial center (groups which have non-trivial elements commuting with all of the other 
group elements), an initial state encoded in one vertex of the network (referred to as reference 
vertex or reference site) can be transferred perfectly to the site labeled by any element of 
the center of the corresponding group. We give explicit analytical formula for the coupling 
strengths in terms of the parameters of the corresponding group association scheme such as 
the diameter d of the scheme and the so-called first and second eigenvalue matrices P and 
Q (which are given in terms of the irreducible characters of the group). By calculating the 
optimal transmission fidelity F opt ., it is shown that for all such networks the perfect transfer 
(Fopt. = 1) can be achieved. We illustrate the method for perfect transfer of a qubit over the 
considered networks with details, and then generalize it to the PST of a qudit over the same 
networks. 

It should be noticed that, the interesting point of this work is that, for the underlying 
graphs of group schemes considered in this paper, the vertices of a the graphs are the elements 
of particular finite groups, where the interactions between qubits associated with vertices are 
governed by the relationship between the vertices defined by the adjacency matrices of the 
scheme; Therefore, according to different hamiltonians (different kind of association schemes 
and consequently different kind of relations or interactions) imposed to a given vertex set, one 
can transfer a given state from a chosen vertex to different vertices. For instance, consider 
a system of 8 qubits each of which has located at a corner of a hypercube; By imposing 
the interactions between these 8 vertices according to the relations of the group scheme over 
G = Z 2 x Z 2 x Z 2 , and taking some coupling strengths between vertices equal to each other, one 



Perfect state transfer 



7 



can obtain a sub scheme which is the same as the Hamming scheme H(2, 3), and transfer state 
of an arbitrary qubit initially prepared at the vertex labeled by |1) = 1 000) to the antipode 
vertex labeled by |8) = 1 1 1 1 ) , perfectly; However, by imposing the other relations defined by 
the other group scheme with 8 vertices, for instance the scheme over the dihedral group D$, one 
can transfer the same state initialized at the vertex labeled by |e) = |000) to the vertex labeled 
by |a 2 ) = |011) in the same graph. In other words, for a given finite set of vertices, one can 
associate different relationships between the vertices by choosing different group association 
schemes in order to transfer a given state from a chosen vertex to different vertices. 

The organization of the paper is as follows: In section 2, we recall some materials about 
association schemes, particularly the group association schemes, their underlying networks and 
the corresponding stratifications. Section 3 is devoted to perfect transfer of a qubit over the 
corresponding networks, where an analytical formula for suitable set of coupling constants in 
particular spin Hamiltonians is given. In section 4, we generalize the method to the perfect 
transfer of a qudit over the same networks. Section 5 contains two examples of the underlying 
networks of group association schemes for which the PST is achieved. The paper is ended with 
a brief conclusion. 

2 Underlying networks of association schemes 

In this section, we review some preliminary tools about the particular networks which are 
considered through the paper. For material not covered in this section, as well as more detailed 
information about association schemes and their underlying networks, refer to [18] and [26]. 
Definition 1 Assume that V and E are vertex and edge sets of a regular resistor network, 
respectively (each edge has a certain conductance). Then, the relations {Ri}o<i<d on V x V 
satisfying the following conditions 
(1) {Ri}o<i<d is a partition of V x V 
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(2) R = {(a, a) : a G V} 

(3) Ri = R\ for < % < d, where R\ = {{(3, a) : {a, (3) G R { } 

(4) For {a, 13) G R k) the number p\- =| {7 G V : (a, 7) G -Rj anrf (7,/?) G | does not 
depend on (a, (3) but only on i,j and k, 

define a symmetric association scheme of class d onV which is denoted by Y = (V, {Ri}o<i<d}- 
Furthermore, if we have = for all k = 0, 2, d, then Y is called commutative. 

For examples of association schemes, consider a cube known as Hamming scheme if (3, 2), 
in which V (the vertex set) is the set of 3-tuples with entries in F 2 = {0, 1}. Two vertices 
are connected if and only if they differ by exactly one entry (see Fig. 1(a)). The distance 
between vertices, i.e. the length of the shortest edge path connecting them, will then indicate 
which relation they are contained in. E.g., if x = (0, 0, 1), y — (0, 1, 1) and z = (1, 0, 1), then 
(x,y) G Ri, (x,z) G Ri and (y,z) G R 2 - 

Let Y = (V, {Ri}o<i<d) be a commutative symmetric association scheme of class d, then 
the matrices A ,A 1 ,...,A d defined by 

1 if (a, (3) G Ri, 



(2-1) 

otherwise (a, f3 G V) 



d 



(A)a,f3 ~ 

are adjacency matrices of Y such that 

A i A j = ^ ii A k . (2-2) 

From (2-2), it is seen that the adjacency matrices A , Ai, ...,A d form a basis for a commutative 
algebra A known as the Bose-Mesner algebra of Y. This algebra has a second basis E ,...,E d 

(known as primitive idempotents of Y) so that 

1 d 
E = —J, EiEj = SijEi, }] Ej — I. (2-3) 

iV i=0 
where, N := \ V\ is the number of vertices (sites) and J is the all-one matrix in A. Let P and 

Q be the matrices relating the two bases for A: 

d 

A i = Y / P ji E j , 0<j<d, 
3=0 
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1 * 



(2-4) 



Then clearly 



PQ = QP = NI. 



(2-5) 



It also follows that 




(2-6) 



which shows that the Pji (resp. Qji) is the j-th eigenvalue (resp. the j-th dual eigenvalue ) of A i 
(resp. Ei) and that the columns of Ej are the corresponding eigenvectors. Thus, m 8 =rank(£ , i ) 
is the multiplicity of the eigenvalue Pji of Ai (provided that Pji ^ Pjk for k ^ i). We see that 
m o — l ? Z)i m i — N, and m 8 =trace£ , i = N(Ei)jj (indeed, E t has only eigenvalues and 1, so 
rank(_E fc ) equals the sum of the eigenvalues). 

Clearly, each non-diagonal (symmetric) relation Ri of an association scheme Y = (V, {Ri} 0<t<d ) 
can be thought of as the network (V, Ri) on V, where we will call it the underlying network of 
association scheme Y. In other words, the underlying network T = (V, R±) of an association 
scheme is an undirected connected network, where the set V and R\ consist of its vertices 
and edges, respectively. Obviously replacing Ri with one of the other relations such as Ri, 
for i ^ 0, 1 will also give us an underlying network T = (V, Ri) (not necessarily a connected 
network) with the same set of vertices but a new set of edges Ri. 

2.1 Stratification 

For an underlying network T, let W = C n (with n — \V\) be the vector space over C consisting 
of column vectors whose coordinates are indexed by vertex set V of T, and whose entries are 
in C. For all f3 e V, let \(3) denotes the element of W with a 1 in the (3 coordinate and in 
all other coordinates. We observe that € V} is an orthonormal basis for W, but in this 

basis, W is reducible and can be reduced to irreducible A-submodules Wi, i = 0, 1, d of the 
Bose-Mesner algebra A (by a A-submodule we mean a subspace Wi of W such that AWi C Wi), 
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i.e., 

W = W ®W 1 ® ...®W d , (2-7) 

where, d is diameter of the corresponding association scheme (for more details see [26]). If we 
define Ti(o) = {(3 E V : (o, (3) E Ri} for an arbitrary chosen vertex o E V (called reference 
vertex), then, the vertex set V can be written as disjoint union of r^a), i.e., 

V=(jT t (a). (2-8) 

i=0 

In fact, the relation (2-8) stratifies the network into a disjoint union of strata (associate classes) 
Tj(o). With each stratum Tj(o) one can associate a unit vector in W (called unit vector 
of i-th stratum) defined by 

V Ki aeTi{o) 

where, \a) denotes the eigenket of ct-th vertex at the associate class Tj(o) and Ki = | (o) | is 
called the i-th valency of the network (/Cj := p% = \{j : (0,7) E Ri}\ = \Ti(o)\). For < i < d, 
the unit vectors of Eq.(2-9) form a basis for irreducible submodule of W with maximal 
dimension denoted by Wo. Since {|0i)}f =o becomes a complete orthonormal basis of Wo, we 
often write 

Wo = X>C|0,>. (2-10) 

i=0 

Let Ai be the adjacency matrix of the underlying network T. From the action of on 
reference state |0o) (|0o) = 1°)) with o E V as reference vertex), we have 

M<h)= E 1/5)- (2-11) 

Then by using (2-9) and (2-11), we obtain 

AM = v^l^>- (2-12) 
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2.1.1 Group association schemes 

Group association schemes are particular association schemes for which the vertex set contains 
elements of a finite group G and the relations Ri are defined by 

R i = {(x,y)\yx- 1 eC i }, i = 0,l,...,d, (2-13) 

where C = {e},C 1 , ...,Cd are the conjugacy classes of G. Then, (G, {Ri}o<i<d) becomes a 
commutative association scheme, and it is called the group association scheme of the finite 
group G. The i-th adjacency matrix Ai is defined as: 

M = Q := £ 9, (2-14) 
where g is considered in the regular representation of the group. Then, we can write 

k=0 

so that the intersection numbers k = 0, 1, d are given by [31] 



jb _ \Ci\\Cj\ ^ x(ai)x(aj)x(ak) , , 

P -~ \G\ \ X (l) {2 - lb) 

where the sum is taken over all the irreducible characters of G. Therefore, the idempotents 

E , ...,Ed of the group association scheme are the projection operators as 

X*(cT> (2-17) 
Thus eigenvalues of adjacency matrices A k and idempotents E k are 



ik 



Qik = d k Xk(ai) (2-18) 



di 

respectively, where dj = XjiX) is the dimension of the irreducible character Xj an d ^k = |Cfe| 
is the /cth valency of the graph. 

It should be noticed that, in the cases that some of the conjugacy classes are not real and 
hence some of the irreducible representations are complex, we encounter with directed under- 
lying graphs or non-symmetric association schemes. In these cases, one can form a symmetric 
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association scheme out of a given non-symmetric association scheme (see the appendix A of 
[26]) so that, the transition probability amplitude between the vertices O £ Co an d /3 G at 
time t is given by 



Wo(*)> 



-i*lXi(«l)« _ 

l^|Z)j^j e di Xi\P) for real representations 

jG| Z)i^i e " ' (Xi(/^) + Xi\P)) for complex representations. 

(2-19) 



for more details, see [26]. 

In the next section, we investigate state transfer over undirected underlying graphs of group 
association schemes, where we will deal with particular hamiltonians which are defined in terms 
of the adjacency matrices of the corresponding association scheme; In the other words, we will 
consider the underlying graphs of group association schemes as spin networks in which the 
interactions betweens spin states associated with the vertices of the graphs are determined by 
the relations Ri of the association scheme defined as in Eq. (2-13) (or equivalently by the 
adjacency matrices A { defined as in (2-14)). 

As regards the above arguments, in the underlying graphs of association schemes, the 
interactions between qubits associated with vertices are governed by the relationship between 
the vertices defined by the adjacency matrices A; L of the scheme; Therefore, according to 
different hamiltonians (different kind of association schemes and consequently different kind 
of relations or interactions) imposed to a given vertex set, one can transfer a given state from 
a chosen vertex to different vertices. For instance, consider a system of 8 qubits each of which 
has located at a corner of a hypercube; By imposing the interactions between these 8 vertices 
according to the relations of the group scheme over G = Z 2 x Z 2 x Z 2 , and taking some coupling 
strengths between vertices equal to each other, one can obtain a sub scheme which is the same 
as the Hamming scheme #(2,3), and transfer state of an arbitrary qubit initially prepared 
at the vertex labeled by ]1) = [000) to the antipode vertex labeled by |8) = [111), perfectly 
(see Fig. 1(a)); However, as we will see in section 5, by imposing the relations defined by the 
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group scheme over the dihedral group D 8 , even by considering only two types of relations (the 
coupling strengths associated with other relations are taken to be zero) (see example 5.1 for 
details), we can transfer the same state initialized at the vertex labeled by |e) = |000) to the 
vertex labeled by |a 2 ) = |011) in the same graph (see Fig.l (b)). In other words, for a given 
finite set of vertices, one can associated different relationships between the vertices by choosing 
different group association schemes in order to transfer a given state from a chosen vertex to 
different vertices. 



3 Perfect state transfer of a qubit over antipodes of un- 
derlying networks of association schemes 

In order to set the scene, let us first recall the algorithm of quantum state transfer over a 
general quantum network, briefly. For more details see Refs. [1, 2, 6, 7]. 

A general finite quantum network is defined to be a simple, connected, finite graph T = 
(V, E), where V denotes the finite set of its vertices and E the set of its edges. A 2-dimensional 
quantum system associated with such a graph is defined by attaching a 2-level (spin-1/2) 
particle to each vertex of the graph so that with each vertex i e V one can associate a Hilbert 
space Hi ~ C 2 . The Hilbert space associated with T is then given by 

H = ® iev Hi = {C 2 )® N , (3-20) 

where N := \V\ denotes the total number of vertices (sites) in T. Now, suppose that we 
impose a particular spin Hamiltonian H (which governs to the interactions between the spin- 
1/2 particles) to the graph so that the total z component of the spin, given by o z tot = J2iev °f * s 
conserved, i.e., [<rl ot ,H] = 0. Hence the Hilbert space Ti decomposes into invariant subspaces, 
each of which is a distinct eigenspace of the operator cr^ ot . Then, for the purpose of quantum 
state transfer, it suffices to restrict our attention to the N- dimensional eigenspace of of ot 
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corresponding to a spin configuration in which all the spins except one are up and one spin is 
down. A basis state for this eigenspace can hence be denoted by the \j), where j is the vertex 
in T at which the spin is down. Thus, {\j) \ j E V} denotes a complete set of orthonormal 
basis vectors spanning the single down subspace. 

The process of transmitting a quantum state from site A to site B proceeds in two steps: 
initialization and evolution. First, a quantum state \iP)a = oi\Q)a + P\1)a £ 7~(<a (with a,f3 E C 
and |a| 2 + |/3| 2 = l)tobe transmitted is created. The state of the entire spin system after this 
step is given by 

\i/>(t = 0)> = |^0...00 B ) = a|0 A 0...00 B ) + /3|1 A 0...00 B ) = a\0) + /3\A), (3-21) 

where, |0) := |CU0...00_b), corresponds to the configuration of all spins up (this state is a zero- 
energy eigenstate of the considered Hamiltonian H). Then, the network couplings are switched 
on and the whole system is allowed to evolve under U (t) = e~ lHt for a fixed time interval, say 
t . The final state becomes 

N 

m o )) = a\0)+PY,fM\j) (3-22) 
j'=i 

where, f jA (t Q ) := {j\e~ iHto \A). Any site B is in a mixed state if |/ab(£o)| < 1> which also 
implies that the state transfer from site A to B is imperfect. In this paper, we will focus only 
on PST. This means that we consider the condition 

\f AB (t )\ = l for some < t < oo (3-23) 

which can be interpreted as the signature of perfect communication (or PST) between A and 
B in time t . The effect of the modulus in (3-23) is that the state at B, after transmission, 
will no longer be but will be of the form 

a|0) + e**/3|l). (3-24) 

The phase factor e l<t> is not a problem because is independent of a and (3 and will thus be 
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a known quantity for the graph, which we can correct for with an appropriate phase gate (for 
more details see for example [1, 2, 15, 16]). 

The model we will consider is an underlying network of a group association scheme over 
a finite group G consisting of N = \G\ sites labeled by {1,2, ...,7V} and diameter d. Then 
we stratify the network with respect to a chosen reference site, say 1 = e (the unit element 
of the group), and assume that the group G has a non trivial conjugacy class, say C m , with 
cardinality one which contains the output site N (i.e., \<f> m ) = \N)). At time t — 0, the qubit 
in the first (input) site of the network is prepared in the state \ipi n ). We wish to transfer the 
state to the Nth. (output) site of the network with unit efficiency after a well-defined period 
of time. The standard basis for an individual qubit (e.g., one can consider state of a spin 1/2 
particle) is chosen to be {|0) = | |), |1) = | |)}, and we shall assume that initially all spins 
point "down" along a prescribed z axis; i.e., the network is in the state |0) = \0a00...00b) ■ 
Then, we consider the dynamics of the system to be governed by the quantum-mechanical 
Hamiltonian 

H = T,Ji E * (<) • ^ (3-25) 
i=o {i,j)eRi 

where, <7j is a vector with familiar Pauli matrices erf, erf and of as its components acting on 
the one-site Hilbert space Hi ~ C 2 associated with the site i e V, and J\ is the coupling 
strength between the reference site 1 and all of the sites belonging to the l-th stratum with 
respect to 1. As it can be easily seen, the Hamiltonian (3-25) commutes with the total spin 
operator (conservation). That is, since we have [cf 0( , H] = 0, hence the total Hilbert space H = 
® iev Hi = (C 2 )® N decomposes into invariant subspaces, each of which is a distinct eigenspace 
of the operator of ot . Then, as it was shown in Ref. [6], one can rewrite the Hamiltonian (3-25) 
in the single down subspace with basis vectors |i) = | f ... f [ | ... |), i — 1, 2, . . . , N as 

i 

follows 

d TV _ 4 d 

H = 2 Ji A + E K i J i L ( 3 - 26 ) 

i=0 ^ i=0 
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Then, by using (2-4), we have 

e -iHt = e -*t V Eto K * J < e" 2i * Eto J ^ E k . 

k=0 

Now, let the m-th stratum of the network contains only one vertex, then one can write 



k=0 

d 



l£ e ^Elo^Q mfc ^E^^^^W- (3-27) 



I I / j lltru I 

M fc=0 \ Lr \ k=0 



In order that the optimal state transfer be achieved, the above probability amplitude must 
take its maximum value. By using (A-i), we obtain 



max 



^U>| 2 = T^ma X \Y l \Qmk\e- it ^ JiP ^ +iekV \ 2 



k 

where, e k is or 1 depending on the sign of Q mk (the sign of Xk{ot m )). The above maximum 
attains if we have 

-tJ2JiPki + e k 7r = (t> + 2l7r -> (PJ) k = £k7V ~f~ 2ln , leZ (3-28) 
i *o 

(0 is a constant phase). Then, by using (2-5), (2-18) and (3-28), the optimal coupling constants 

(for which the optimal state transfer is achieved) are obtained as 

J k = T77]Z^ 7 Qki = j^T 2^ 7 diXi{a k ). (3-29) 

so that the optimal fidelity is given by 

cZ 1 ^ 

F opt . = max(0 m |e" im |0 o ) = ^HlQmil = 777, £<*i|Xi(om)l- O 30 ) 

M i=o M i=0 

From the fact that a m G C m belongs to the center of the group (and so commute with all 
elements of the group), the Schur's lemma implies that the irreducible representation of a m 
is equal to al for some ji e C. Now, let a be the order of a m E G {a a m = 1). Then, clearly 
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fj, must be the ath root of unity, i.e., we have /x = e 2m l a and so = 1. Then, we will have 
IXfc^m)! — \n\dk — dk and the optimal fidelity attains 

Fo P t=^jZdl = l. (3-31) 

M k=0 

where we have used the well known identity J2k=o d\ = \G\ from group theory. The above result 
indicates that all finite groups with non-trivial center can allow state transfer with optimal 
fidelity one, i.e., perfect state transfer. 



4 Generalization to perfect state transfer of a D-level 
quantum state 

In the D N dimensional Hilbert space associated with system of N, D-level quantum states 
(H ~ (C D )® N ), the Hamiltonian (3-25) can be generalized as 

H = jlJi E (4-32) 

i=o (i,j)eR, 

where, A, is a D 2 — 1 dimensional vector with generators of SU (D) as its components acting 
on the one-site Hilbert space Hi ~ C D . 

Let us denote a state in which the i-th site has been exited to the level v by |z/j) = 
|0 0^z^0...0). Then, the hamiltonian H can be diagonalized in each subspace 

i 

spanned by the vectors |^), % = 1, . . . , N, for v = 1, D — 1 (for more details see [7]). If we 
call the states with only one site excited, as one particle states and the subspace spanned by 
these vectors comprise the one particle sector of the full Hilbert space, then, the whole one 
particle subspace S can be written as 

In the other words, in D N dimensional Hilbert space H, we deal with D — 1 one particle 
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subspaces (recall that, each of these subspaces has dimension N). In the case of a system of 
N qubit (D — 2), we have only one one-particle subspace of dimension N. 

As it was shown in Refs. [6], [7], the hamiltonian (4-32) restricted to each one particle 
subspace S^ u \ for v — 1, 2, . . . , d — 1, can be written in terms of the adjacency matrices A\ as 

H = 2 £ JA + (N - 2) ° ~ AN £ J,«,l. (4-33) 

Then the same procedure can be applied in order to obtain the optimal fidelity. The Eq.(4-33) 
is similar to the Eq. (3-26) apart from the non-important second term which is multiple of 
the identity. Therefore, maximizing the probability amplitude (0 m |e~* m |0o) gives the same 
constraints on coupling constants J\ and leads to the same optimal fidelity (3-31). 



5 Examples 

In this section, we consider two examples of group association schemes and investigate PST 
over them. 

5.1 The dihedral group D<m 

The dihedral group G = D 2n is generated by two generators a and b with the following relations: 

D 2n = (a, b : a n = 1, b 2 = 1, b^ab = a' 1 } (5-34) 

We consider the case of even n = 2m, the case of odd n can be considered similarly. The 
Dihedral group D 2n with even n = 2m has m + 3 conjugacy class so that C m = {a m }. For 
even n = 2m, the m + 3 conjugacy classes are given by 

C = {1}, C r = {a r ,a- r } ; 1 < r < m-1, C m = {a m }, C m+1 = {a 2 >b ; j = 0, . . . , m-1}, 

C m+ 2 = {a 23+1 b ; j = 0,...,m-l}. 
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The character table is given by 



Dor, 


e 


a m 


a r (1 < r < ttt. — 1) 


b 


ah 


AO 


i 


i 


i 


i 


i 

± 


Xi 


1 


1 


1 


-1 


-1 


X2 


1 


(-l) m 


("l) r 


1 


-1 


X3 


1 


(-l) m 


(-iy 


-1 


1 


Ipj (1<J<?7Z — 1) 


2 


2(-iy 


2 cos(2ir jr /n) 









Then, by using the result (3-30) we obtain 

Y rn+2 -y 

Kpt. = Tq\zZ d k \Xk(a m )\ = — {4 + 4(m- 1)} = I. 

I I fc=0 

From (3-29), one can evaluate the optimal coupling constants J\. For example, for n = 4 we 
obtain 

7T 2"7T 

Jo = J2 = J4 = 0, Ji = — , J3 = — , 

zt to 

where the above values are obtained by choosing = 7r, l — 3, l\ — — 1, Z 2 = —^3 = 2 and 
/ 4 = in the Eq. (3-29). It is seen that by choosing = vr, all coupling constants become zero 
except for J x and J 3 (see Fig. 1(b)). 



5.2 The Cubic group T h 

One of the point groups with high or polyhedral symmetry is the cubic group Th of order 24 
with pyritohedral symmetry. This group is isomorphic to At x C 2 , where A A is the alternating 
group of order 12. It is the symmetry of a cube with on each face a line segment dividing 
the face into two equal rectangles, such that the line segments of adjacent faces do not meet 
at the edge. The symmetries correspond to the even permutations of the body diagonals and 
the same combined with inversion. The group Th has 8 conjugacy classes with the following 
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character table 



T h 


C 


Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


CV 


Xo 


1 


1 


1 


1 


1 


1 


1 


1 


Xi 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


X2 


1 


1 


1 


1 


IV 


tv 2 


IV 


u 2 


X3 


1 


1 


1 


1 




tv 


v 2 


V 


X4 


1 


-1 


-1 


-1 


tv 


iV 2 


—v 


-v 2 


Xb 


1 


1 


-1 


-1 


CO 2 


tv 


-v 2 


—v 


X6 


3 


-1 


3 


-1 














X7 


3 


-1 


-3 


1 















with tv := e 27 ™/ 3 . By combining the classes C 4 and C 5 and denoting it by C4 = C4UC5, and C e 
with C7 and denoting the new obtained class as C 5 , one can obtain a symmetric association 
scheme with undirected underlying graph. Now, using the result (3-29) one can evaluate the 
optimal coupling constants J\ for, / = 0, 1, . . . , 5 as follows: 

7T 

Jo — J\ — J2 — J4 — 0, J3 = J5 = 7—, 

where, we have chosen = 7r/2, l = —l\ = 3 and l 2 = h = h = h = in the Eq. (3-29). By 
these choices, all coupling constants become zero except for J 3 and J 5 (see Fig. 2). 

5.3 The Clifford group 

The Clifford algebra with n generator matrices 71,72, •••,7n, obeys the following relations [32] 

lilj + Ijli = (5-35) 

Thus, the 7's have square 1 and anti-commute. The Clifford group denoted by CL(n) has 2 n+1 
elements as 

CL(n) = {±l,±7 il ...7j.; % x < ... < = 
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where, i r G {1,2, ... ,n}. We suppose n > 2 throughout. It is well known that [32], the 
center of CL(n) denoted by Z(CL(n)), consists of {±1} if n is even and {±1, ±7!...7„} if n 
is odd. CL(n) has 2 n one-dimensional representations, each real. In each such representation, 
U(—l) = I; Any irreducible representation with dimension greater than 1 has U(— 1) = — /. 
For even n, the conjugacy classes are given by 

C = {1}, Ci = {-1}, C 2 = {71, -7i}, • • • , Cj = {7u-7i,-, -7u-7iJ, 

C 2 n +1 = {71. ..7„, -7i-7i}> 

whereas for odd n, we have 

C = {1}, Ci = {-1}, C 2 = {71, -7i}, • • • , Cj = {7u-7i,-, -7<i-7<j}, 

CWl = {7l —Tn}, ,C2™+2 = {— 7i—7i}- 

In the following we consider only the case of even n, the case of odd n can be considered 
similarly. 

The characters of the 2 n one dimensional representations are given by 

Xfc(l) = Xfc(-l) = 1, 0<fc<2 n -l, 
X2n {± lA ) = ±5 A% 2 n l 2 X2n {l) = 2 n l\ X2 „(-l) = -2"/ 2 . 
Then, by using the result (3-30) we obtain 

F op t. = 2^rEd*lx fc («m)l = ^ T {2" + 2"/ 2 (2«/ 2 )} = 1. 

Even though we considered only two examples of PST over underlying networks of group 
association schemes in details, for all underlying networks of group association schemes for 
which the corresponding group has non-trivial center, PST can be achieved. For example, 
underlying networks of group association schemes associated with abelian finite groups and 
all p-groups (a p-group is a group whose order is a power of the prime number p; the centers 
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corresponding to the p-groups are non-trivial [31]) allow PST. It should also be noticed that 
for direct product groups G = G 1 x G 2 , so that PST can be achieved for Gi and G 2 at the same 
time t , PST will be achievable for the group G, since the transition amplitudes (0fc|e~ li? *|0o) 
corresponding to the underlying network associated with group association scheme over G 
are product of the amplitudes corresponding to the groups Gi and G 2 (for more details see 
[26]). Particularly, for any finite group G with non-trivial center, not only PST over un- 
derlying network of group association scheme over G is achieved, but also it is achievable for 
direct product group G x G x . . . x G . For instance, the n-dimensional hypercube network (n- 

n 

cube) can be viewed as the underlying network of group association scheme over the product 
group Z 2 x Z 2 x . . . x Z 2 and PST over it can be achieved by only nearest neighbor coupling 

n 

strengths. 

6 Conclusion 

Perfect state transfer of a qudit over antipodes of the so called underlying networks of group 
association schemes, was investigated, where by using the group properties and the algebraic 
structure of these networks (such as Bose-Mesner algebra), an explicit analytical formula for 
coupling constants in the Hamiltonians (in terms of the irreducible characters of the corre- 
sponding group) was given, so that the state of a particular qudit initially encoded on one site 
can be perfectly evolved to the opposite site without any dynamical control. 
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Appendix A 

As regards the arguments of section 3, for the purpose of the perfect state transfer, we con- 
sider the underlying graphs of group association schemes in which the corresponding group 
has non-trivial center in order to the group have no-trivial one element conjugacy classes or 
stratas with size one, i.e., the last stratum of the graph denoted by T m (o) contains only one 
site (K m = |r m (o)| = 1). Then, by imposing the constraint that the transition probability 
amplitude (0 m |e~* m |0 o ) be equal to an arbitrary phase such as e l6 and consequently the tran- 
sition probability amplitudes {4>i\e~' lHt \(f) ) be zero, for all % ^ m, one can achieve PST from 
the reference site |0o) = |e) associated with the first stratum to the site of the last stratum 

\<t>m)- 

Even though, as it was shown in section 3, we need to evaluate only the probability am- 
plitude ((p m \e~ tHt \(po) in order for the corresponding fidelity to be maximized and PST be 
attained, here we give the transition probability amplitudes between the reference site |e) and 
any group element (3 G G, in terms of the coupling constants and irreducible characters of the 
group which can be used in some areas of quantum information such as generating quantum 
entanglements between any two sites of a given finite network. To do so, given a finite group 
G, let (3 e G belong to the l-th conjugacy class of the group. In general, assume that the 
group G has R real conjugacy classes (Cj = C" 1 , for % — 0, 1, . . . , R— 1) and so it will posses R 
real representations; Now, we use the stratification and algebraic properties of the underlying 
graph of the group association scheme over the group G, to write 

(/3|e-^|0 o > = -U(^|e-^U> = ~^{j:e- ite ^d k xM+ E e^d*(Xfc(ai)+Xfc(ai))}, 

V K l \ Lt \V k I k=0 k=R 

(A-i) 

where 

R 1 d+l-R 

© Xfc := E + E ^(XM + xMi)) (A-ii) 

1=0 a k l=R a k 

where we have used the fact that Xk(®i) = Xk( a i)- As it is seen from (A-ii), we have 6 Xfe = 6 x - fe . 
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It should be also noticed that, for groups G with non-real conjugacy classes, the corre- 
sponding group association schemes are not symmetric and so the corresponding underlying 
graphs are not undirected; Therefore, for the purpose of state transfer, we need to symmetrize 
the group schemes in order to obtain undirected graphs. Now, we recall that for any non real 
conjugacy class (Cj ^ Cf 1 ), Cf 1 is also a conjugacy class of G with the same size of Cj, i.e., 
\d\ = IC" 1 !; Then, as in Ref. [26] was shown, the symmetrized classes C,, for % — 0, . . . , R — 1; 
Ci = Ci U C" 1 , for i = R, . . . , d+1 ~ R , form a symmetric association scheme which is a sub 
scheme of the non-symmetric group association scheme over G. Now, by choosing the identity 
element of the group as reference vertex and stratify the underlynig graph with respect to 
it, the resulted symmetrized sub scheme can be considered for the purpose of PST if its last 
stratum denoted by V m (e) has size one, i.e., we must have \C m = C^l = 1 in order for PST 
to be studied. 
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Figure Captions 

Figure 1: (a) Denotes a vertex set {1, 2, ... , 8} with interactions according to the relations 
of Hamming scheme H(3, 2), where with only one non-zero coupling strengths J±, PST between 
the vertex (000) and the vertex (111) is achieved, (b) Denotes the same vertex set {1,2,..., 8} 
with interactions according to the relations of the group association scheme over dihedral group 
Dg, where with only two non-zero coupling strengths J± and J3, PST from the vertex e to the 
vertex a 2 is achieved; The solid lines denote the interaction coupling J x and the dashed lines 
denote the coupling strength J 3 . 

Figure 2: Shows the underling graph of the symmetrized scheme obtained from the group 
scheme over the Cubic group T h , where non-zero coupling strengths J 3 = J 5 and J = J\ = 
J 2 = J4 = 0, PST from C to C 2 is achieved. 



